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Abstract 

We present a comprehensive description of vector chromatography that includes deterministic 
and stochastic transport in ID-periodic free-energy landscapes, with both energetic and entropic 
contributions, and highlights the parameters governing the deflection angle, i.e. the Peclet number 
and the partition ratio. We also investigate the dependence of the deflection angle on the shape 
of the free-energy landscape by varying the width of the linear transitions in an otherwise di- 
chotomous potential. Finally, we present experimental results obtained in a microfluidic system in 
which gravity drives the suspended particles and, in combination with a bottom surface patterned 
with shallow rectangular grooves, creates a periodic landscape of (potential) energy barriers. The 
experiments validate the model and demonstrate that a simple, passive microdevice can lead to 
vector separation of colloidal particles based on both size and density. 



Microfluidic systems for chemical and biological separation have shown great promise and 
opened the door for exciting new technologies. A number of separation systems based on 
driving suspended particles through a periodic stationary phase, for example, take advantage 
of the unprecedented control on the geometry and chemistry provided by micro-fabrication 
techniques. In fact, one- dimensional (ID) separation along the driving direction has been 
demonstrated in different microdevices with periodic stationary phase, ranging from en- 



tropic trap arrays [1] to asymmetric structures acting as a ratchet |2j. The description of 
ID transport of suspended particles past periodic entropy barriers, and to a lesser extent 
energy barriers, has also received considerable attention and rigorous results are available 



for the effective mobility of single particles 



13( | . Two-dimensional (2D) separation meth- 



ods, in which the different species constituting a sample migrate in different directions, 
enabling their continuous fractionation and, in general, providing greater selectivity than 
ID techniques, have also been developed based on periodic stationary media, and have been 
categorized as vector chromatography (VC) jl4j ]. Notably, VC can be obtained in planar 
microfluidic devices via a straightforward extension of the aformentioned ID methods, by 
driving the particles at an oblique angle with respect to the periodic direction. This is the 
case when suspended particles are driven through force fields that are periodic in one of 



the directions of the separation plane and invariant in the other [l5KL7j|. Although a case- 
by-case analysis in the deterministic limit provided good agreement with experiments in 
these systems, a general description is lacking. In this letter, we extend previous results to 
obtain a comprehensive description of planar VC in terms of the ID-periodic free-energy of 
the system, including energetic and entropic contributions, that captures the deterministic 
and Brownian limits. This unified description highlights the key parameters governing the 
migration angle of different species and their relevance to the design and optimization of 
fractionation devices. We also present experimental results obtained in a microfluidic sys- 
tem in which gravity not only drives the particles but also, in combination with a bottom 
surface patterned with shallow rectangular grooves, creates a periodic landscape of potential 
energy barriers. The experiments agree well with the theoretical results, exhibit several of 
the qualitative features predicted by the model, and show the separation capability of the 
device. 

Consider the transport of non-interacting Brownian particles in a planar microfluidic 
device, driven by a constant external force F (oriented at an angle Op = aictan(F y /F x )) and 



moving through a potential energy landscape V(x, y) that is periodic in one of the directions 
of the separation plane, say x axis [periodic direction), and invariant along the other, say y 
axis (uniform direction). Let us assume that the driving force is contained in the separation 
plane, with any vertical component conveniently incorporated into the potential V. The 
asymptotic distribution of particles in a unit cell is given by the steady-state solution of the 



Smoluchowski equation for the reduced probability density P°°(x) [18l. Il9j. 



= V • J(x) = V • ( U(x)P°°(x) - D(x) ■ WP°°(x)) , (1) 

where J(x) is the probability density flux, U(x) is the instantaneous particle velocity and 
D(x) is the diffusion tensor. In the low Reynolds number limit the velocity of the particle is 
a linear combination of the forces acting on it, U(x) = M(x) ■ (F — W(x)), where M(x) is 
the mobility tensor, which locally satisfies the Stokes-Einstein relation D(x) = ksTM(x). 
P°°(x) is periodic in x, satisfies the no-flux condition in z, and is normalized, J T P°°(x)dV = 
1, where r is the volume of the unit cell. Given P°°(x) it is straightforward to compute 
the components of the average migration velocity U XiV = J T J x>y dr and the migration angle 
9 = arctan([7 y / 'U x ) , which is the relevant parameter in VC, by means of macrotransport 
theory 



In planar microfluidic devices the particles are usually highly confined in the vertical 
direction, either geometrically, as in the case of entropic trapping, or due to particle-wall 
interaction potentials with narrow secondary minima {?], [ill 13- m ^ ms case > it is valid to as- 
sume fast equilibrium in the cross-section, which in this context is known as the Fick- Jacobs 
(F J) approximation j^, I201 2 



fast degrees of freedom 



]23 



22] and is analogous to other projection methods that eliminate 
3-25]. It is then possible to write the probability distribution in 
terms of the marginal probability density, V(x), and the equilibrium conditional distribution 
in the cross-section, p(z\x) eq , 

P°°(x) w V(x)p(z\x) eq = V(x)Q- 1 e-^ x ' z \ (2) 

where (3 = (ksT)^ 1 and Q(x) = J exp(—(3V(x,z))dzdy is the local partition function. In 
this approximation, the average migration velocity in the uniform direction is determined by 
the average mobility, that is U y = (^Jq x {M yy ) eq V(x)dx^j F y , where M yy is the hydrodynamic 
mobility and, for any function f(x, z), (f) eq = I f(x, z)p(z\x) eq dzdy is the local equilibrium 
average over the cross-section. Before we calculate the average velocity in the periodic 



direction, we note that U x = Jq x J x dx = £ X J X , where J x = J J J x dydz is the total flux 
through any cross-section in the periodic direction and is constant in steady state. Then, 
integrating Eq. (OQ) over the cross-section and using the FJ approximation we obtain 

J* = (M xx ) eq { IF* ~ ?{x)] V - k B T^ , (3) 

where J-(x) is the mean force due to the potential, and following Zwanzig's derivation for 
the purely diffusive case (3]], we write it in terms of the local free-energy of the system, 
A(x) = —ksT In Q(x), 

It is clear from Eq. fl3]) that the total flux in the periodic direction, and therefore U x , 
have both diffusive and convective contributions. The first convective term shows that, 
as expected, even in the absence of a periodic potential (i. e. J 7 = 0), a macroscopic 
anisotropy in the mobilities, {M xx ) eg 7^ (M yy ) eq , could lead to a non-zero deflection angle, 
A9 = 9 — 9p 7^ 0. Furthermore, the dependence of the average mobility on x can also result 
in a non-zero deflection angle, independent of the local isotropy of the mobility tensor, given 
that the contribution of the diffusive flux to U x (last term in Eq. fl3])) would not vanish 
in this case [13]. For simplicity, however, we shall assume that the mobility functions are 
constant and equal, M xx = M yy = M. In this case, we have U y = MF y and it is clear 
that a non-zero deflection angle implies that U x 7^ MF X , and vice versa. Moreover, in this 
case the diffusive contribution to U x vanishes, and therefore only J-{x) can contribute to a 



non-zero deflection angle. Solving Eq.([3]) [13|, |23] we obtain 

Pe fi _a MS * P+ 1 



tan 9 = tan 9 F 



l-e- pe 7o 



■[dxe-WW jdie^ , (5) 

JO Jx i 



where (3A'(x) = (3A(x) — Pex, and x = x/£ x . The Peclet number, Pe = (3F X £ X , measures 
the relative magnitude of convective to thermal transport in the system and is one of the 
two dimensionless parameters dictating the behavior of the migration angle. Alternatively, 
one can use a characteristic value of the mean force, such as its maximum value J- max , and 
consider the normalized driving force / = F x /J r max as an independent parameter. However, 
we shall see below that the Peclet number and the normalized force are complementary for 
the description of a given system, in that one is the appropriate parameter to consider when 
the other diverges. The second parameter is the partition ratio /C, which measures the spatial 



variations in the distribution of particles in equilibrium. Its importance becomes apparent 
when we adopt the amplitude of variations in the free-energy, AA, as the characteristic 
energy scale to write (3A(x) = (In K) A{x) , with K = exp(/3ZL4), which in the context of 



26). 



transition-state theory, corresponds to an Arrhenius factor 

In order to investigate the role of these parameters on the migration angle we consider 
a cosine potential, A{x) = l/2cos27nc, and a dichotomous potential with linear transitions 
(LTD potential), given by regions of constant potential, A = and A = 1, connected by 
linear transitions of width 5 (see inset in Fig. [2]). The LTD potential in the limit 5 = 
corresponds to a square wave (SW) potential. In all cases, the effect of the periodic 
potential is to reduce the average velocity in the periodic direction j^, U x , resulting in 
positive deflection angles (for the cases shown here in which Op = 45°, this implies that 
0° < AO < 45°). Fig. [1] shows the deflection angle as a function of Pe for the SW (solid 
curves), LTD (dashed lines), and cosine (dotted line) potentials. It is clear that the deflection 
angle decreases with Pe and increases with /C (arrow direction), independent of its entropic 
or energetic origin. Fig. [1] also shows that for a given K, and Pe, the smaller the transition 
region in the LTD potential the higher the deflection angle, with, as expected, A9 converging 
to the curve for the SW potential in the limit 5 — > (arrow direction). In Fig. [2] we plot the 
deflection angle as a function of the normalized force and consider the effect of Brownian 
motion for the different potentials. Specifically, we compare the deflection angle obtained 
at a finite partition ratio for the cosine and LTD potentials with that in the deterministic 
limit (Pe — > oo and finite /). Note that Figs. [1] and [2] are complementary, in that they allow 
us to investigate independent limits, i. e., the SW potential limit (/ — > and finite Pe) and 
the deterministic limit (Pe — > oo and finite /), respectively. In all cases, the deflection angle 
decreases with /, which is similar to the behavior observed with Pe in Fig. [TJ Fig. [2] also 
shows that Brownian motion allows the particles to cross the potential barriers even when 
/ < 1. This is in contrast to the deterministic case, in which particles are locked to move 
along the uniform direction. We also investigate the effect that the transition region S has 
on the deflection angle. In the deterministic limit, the deflection angle for the LTD potential 



has the simple analytical expression 

tan# 



28| 



1 f6 f6 



(6) 



tan Op 

and it is clear that larger transition regions lead to larger deflection angles for any / > 1 

5 




l-4tanh(Pe/4)/Pe 



FIG. 1. Deflection angle as a function of the Peclet number. Solid lines correspond to the SW 
potential with log/C = f3AA = 1,2,3,4 and 6. The dashed lines correspond to the LTD potential 
with (3AA = 3, transition regions 5 = 0.01,0.1,0.2,0.3,0.5, and e\ = ei = 0.5 — 5. The arrow 
traverses curves of increasing /C for the SW potential and curves of decreasing 5 for the LTD 
potential. The evolution (a — > b) of a purely entropic system upon a temperature increase is 
shown. 

The reason is that as 5 increases (at constant /), the particle is deflected for a longer time 
as it crosses the potential barrier. Note the difference in the contributions coming from 
the regions with ±/ to the transit time, as shown by the corresponding terms in Eq. [61 

6 



In the presence of Brownian motion we observe the same trend for large driving forces 
/, as expected. On the other hand, as the driving force decreases and barrier hoping is 
dominated by thermal motion, the behavior reverses and larger transition regions lead to 
smaller deflection angles. This crossover between the deterministic and Brownian cases as 
/ decreases is consistent with the behavior observed as a function of Pe in Fig. [TJ In fact, 
the limits Pe <C 1 and / < 1 coincide in the linear response regime, where the reduction 
in mobility (and diffusivity) is given by (/ Qdx J Q~ 1 dx) [29| . Fig. [3] shows the effect of the 
partition ratio on the deflection angle for the LTD potential with a given 5. Clearly, the 
deflection angle increases as the partition ratio increases, converging to an asymptotic curve 
for /C — > oo. This upper limit coincides with the deterministic limit for a purely energetic 



potential of mean force j9|. In terms of separation devices, it becomes clear that in order to 
obtain large deflection angles and high selectivity it is desirable to operate around / ^ 1. 

The results presented in Figs. HH3] also reveal the role of temperature in different sep- 
aration systems ({5|. In entropic trapping (AA oc kT), for example, the partition ratio is 
determined by the ratio between the available configurations in the slit and well regions js], 
and is thus independent of temperature. On the other hand, both Pe and / decrease with 
temperature, thus leading to larger deflection angles. In fact, increasing the temperature 
for a given geometry of the entropic traps corresponds to the system moving along curves 
of constant partition ratio, as shown in Figs. [1] and [3] On the other hand, in the purely 
energetic case AA and / are independent of temperature. Therefore, for a given potential 
energy landscape, increasing the temperature reduces the partition ratio and the deflection 
angle. This corresponds to the system moving along curves of constant / and decreasing 
partition ratio, that is vertical lines in Fig. [31 as shown. 

We performed experiments in a microfluidic system in which suspended particles are 
driven over a glass surface on which we had (isotropically) etched rectangular grooves in 
parallel. The periodicity is 20 /zm, the shallow grooves are A7-L = 65 nm deep and 13 fim 
wide. The channel is high enough to neglect confinement effects. Gravity induces periodic 
energy barriers due to the presence of grooves, and also drives the particles by tilting the 
bottom surface an angle Qf We use silica particles (4.32 /im and 2.14 //m diameter) and 
polystyrene particles (4.31 \xm diameter). The gravity-induced partition ratio for a particle of 
radius a is given by K = exp (4/37ra 3 Z\p g cos9 t Al-L/k B T) where Ap is the buoyant density 



of the particles and g is the acceleration due to gravity 



13|. For 6 t = 0° the partition 




f=F/T 



FIG. 2. Effect of 5 on the deffection angle as function of the normalized force in the stochastic and 
deterministic regimes. The solid lines correspond to the cosine potential while the dashed lines 
correspond to the LTD potential for different transition regions and e\ = €2 = 0.5 — 5. The curves 
corresponding to the Brownian case and the respective deterministic curves have the same line 
style. The inset shows schematics of the cosine and the LTD potentials. 

coefficient of the 4.32 fim silica particles is at least two orders of magnitude larger than that 
for the 2.14 /zm silica and 4.31 /zm polystyrene particles. Thus, at small tilt angles, the 4.32 
fim silica particles should experience much larger deflections than the other particles, which 
would demonstrate that it is possible to fractionate particles by size or density. In Fig. H] 
we show the experimental results for the deflection angle,for a forcing angle dp = 45°, as a 
function of the tilt angle (note that unlike Pe and /, the tilt angle is common to all particles 
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f=F/T 



FIG. 3. Effect of the partition coefficient on the deflection angle as a function of the normalized 
force for the LTD potential for a given width of the transition region and e± = 62 = 0.5 — 5. 
The arrow traverses curves of increasing partition ratio. The dashed curve corresponds to the 
deterministic limit. Upon increasing the temperature, a system at point a "moves" along the same 
curve toward point b for a purely entropic potential, and vertically down toward point c for a purely 
energetic potential. 

in a given experiment). The theoretical curve corresponds to the LTD potential, with S 
calculated from the best fit to the experimental data and representing an effective transition 
region in the interaction between a suspended particle and the bottom grooves. We obtain 
good agreement for the 4.32 /im silica particles with 5 — 0.10 ± 0.01 (2.0 ± 0.2 /zm), which 
compares well with an order of magnitude estimate S zoi = 2y / 2(a + h e ) {AH + k^ 1 ) = 1.4 

9 




FIG. 4. Deflection angle as a function of the tilt angle. The circular (dashed) and triangular 
(dotted) symbols (curves) correspond to the experimental results (SW potential) for the 4.31 um 
polystyrene particles and to the 2.14 um silica particles, respectively. The solid dark line is a fit of 
the experimental data to the results of the LTD potential for the 4.32 /jm silica particles (square 
symbols) using the width of the transition region as a fitting parameter (5 = 0.10), while the light 
solid line corresponds to the respective deterministic limit. See supplementary information for the 
standard deviations and for a video showing a representative experiment for the 4.32 um silica 
particles. 



lim. obtained by extending the concept of a zone of influence [30] for a particle suspended at 
its equilibrium separation from the wall (h e = 259 nm) and in the vicinity of a step 28]. The 



theoretical curves for 2.14 /im silica and 4.31 /im polystyrene particles are insensitive to the 
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width of the transition region, with differences smaller than 2.5°, and the data is compared 
to SW potentials, with good agreement. In these latter cases, the particles easily overcome 
the energy barriers due to thermal fluctuations, significantly reducing confinement effects 
and leading to small deflection angles. These experiments show that gravity, along with a 
bottom surface patterned with slanted periodic grooves, could be used to separate particles 
according to their mass. More generally, the role of gravity can be effected or enhanced by 
other fields, e.g. electric, dielectrophoretic, magnetic, etc., potentially leading to a versatile 
separation technique. 
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